Dilatation symmetry of the Fokker-Planck 
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^ Based on the canonical formalism, the dilatation symmetry is implemented to the 

a 

q Fokker-Planck equation for the Wigner distribution that describes atomic motion in an 

> 

• i-H 

^ optical lattice. This reveals the symmetry principle underlying the recent result obtained 

u ' " 

by Lutz [Phys. Rev. A 67, 051402(R) (2003)] on the connection between anomalous 
transport in the optical lattice and Tsallis statistics in the high-energy regime. Lutz's 
discussion is generalized to the nonstationary case, and the condition, under which the 
solution distribution decays as a power law, is derived. 
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Anomalous diffusion [1,2] has continuously been attracting attention over decades. It 
can be observed, for example, in turbulent flows [3], dissolved micelles [4], chaotic 
dynamics [5], porous glasses [6], and subrecoil laser cooling [7]. 

Anomalous diffusion, or superdiffusion more specifically, may be a signal of scale- 
free non-Gaussian statistics. Levy statistics [8] is one such example and is characterized 
by a stable class of distributions. In the one-dimensional symmetric case, the Levy 
distribution indexed by y may be defined through the characteristic function of the 

stretched-exponential form: L y (x) = (1 / 2k) ^dk exp (ikx - a I k 1 7 ) , where y e (0, 2) 

and a is a positive constant. Excepting the Gaussian case ( y — > 2 - 0), L y (x) decays as 
a power law: L y (x) ~l x I" 1 " 7 . An important point is the following. There exists a 
mathematical result termed the Levy-Gnedenko generalized central limit theorem [9], 
which states that, by /V-fold convolution, a distribution with divergent lower moments 

tends to one of the Levy stable class, \h y (*)} , in the limit /V — > °° if such a limit is 

convergent. This is in contrast to the ordinary central limit theorem for distributions 
with the finite second moments in normal Gaussian statistics. Now, assume L y (x) to 
describe a single jump in the probabilistic process. After /V independent jumps, one has 
L i y N) (x) = N~ Vy L y (xl N X ' Y ). Identifying N with time, t , one obtains the scaling of the 
distribution, P(x, t) = L ( y N> (x) = t~ 1/y P(x I t Uy ), implying that spreading of the 
distribution follows the law of superdiffusion, ~ t 1 ' 7 , which should be compared with 
the case of normal diffusion, ~t U2 . 
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In a recent paper [10], Lutz has reexamined the problem of anomalous transport of 
atoms in an optical lattice. He has considered the following generalized Fokker-Planck 
equation for the marginal Wigner distribution of the momentum, p, of an atom in the 



optical potential: 



at dp dp 



D(p) 



dWjp, t) 
dp 



(1) 



Both the drift, K, and the diffusion coefficient, D, explicitly depend on the momentum 



of the atom as follows: 



K(p) = - 



ap 



i+(p/p c y 



D(p) = D 0+ ° l 

l + (p/ p c ) 



(2) 



Here, a and p c are the damping coefficient and the capture momentum, respectively. 
D Q and D { are constants related to fluctuations caused by the photon processes. Eq. (1) 
can be derived, after spatial average, from the master equation for the full Wigner 
distribution constructed by quantum dynamics of the atom-laser interaction (see Ref. 
[10] and the references therein). It has been noticed [10] that K and D satisfy the 



relation 
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= ^-lne q (-Ps(p)), 
a p 



(3) 



where e q (s) is the "^-exponential function" defined by e q (s) = (1 + (1 - q)s) l ^ (l q) with 
the notation (a) + = max {0, a} , and 



2£> n „ a 



9 = 1 + —^-, P = — — , e(p) = p 2 . (4) 

ap c 2(D + D,) 



In the limit g— >1, e q (s) converges to the ordinary exponential function, e'.FromEq. 
(3), it follows that the exact stationary solution of Eq. (1) is given by [10] 



W q (p) = Y^e q (-pe(p)) (5) 



with 



Z q (f3) = jdpe q {-(3e(p)), (6) 



which optimizes the Tsallis entropy [11,12], S q [W] = (1 - q)~ l [J dp W"(p)-lj, under 

the constraints on normalization of W and the "^-expectation value" of the "energy", 
e(p) [13,14]. [3 is related to the Lagrange multiplier associated with the ^-expectation 



value of e(p) (see Refs. [13,14] for more details). From these discussions, we see that 
the entropic index, q, is determined by dynamics as in Eq. (4). 

It should be mentioned that since q is larger than unity, the distribution in Eq. (5) 
decays as a power law. In fact, it has explicitly been shown [15] that if q e (5 / 3, 3), 
many-time convolution of the distribution of the form in Eq. (5) converges to the Levy 
distribution with the index, y - (3 - q) I (q - 1), in accordance with the Levy-Gnedenko 
generalized central limit theorem. 

Now, emergence of an asymptotically power-law distribution from the linear Fokker- 
Planck equation is rather peculiar [16-18]. In the present case, the origin of the power 
law is in the behavior of the ratio in Eq. (3), that is, 



(7) 

Dip) p 



in the high-energy regime. One may wonder if there is an underlying principle for the 
emergence of this scale-free nature. In what follows, we reveal such a principle by 
taking advantage of the dilatation symmetry implemented to the linear Fokker-Planck 
equation. As a result, without assuming stationarity, we shall obtain the condition for 
the solution, W, to be asymptotically scale-free, which turns out to be more general than 
Eq. (7). 

Our starting point is the variational principle for the kinetic equations [19-21]. The 



action and the Lagrangian density respectively read 



I[W, A] = \\dt dp £ (W, A, dW I dt, dAI dt,dW I dp,dAI dp), (8) 



£-1 

2 



dt dt J dt dp dp' 



where A(p, t) is an auxiliary field. Performing integration by parts and dropping the 
boundary terms, we see that Eq. (8) can also be expressed in the following form: / = 

-jdt (dA/dt + KdA/dp + d(DdA/dp) /dp), where (A) stands for the ordinary 

expectation value of A: (A) = | dp A(p, t) W(p, t). Taking the variations of the action 

with respect to A and W, we obtain Eq. (1) and dA/dt + KdAldp 
+d(DdA/dp)/dp = 0, respectively. 

Let us proceed to developing the canonical formalism. The canonical momenta 
conjugate to W and A are given by 



d£ 1 

d(dW/dt) ~ 2 



n ff = ^,„„„ =7A, (io) 



d£ 1 

d(dA/dt)~ 2 



n A = ^ t — = --w, (ii) 
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respectively, leading to a pair of the constraints 



*i = n w --A*o, 



x 2 = n A + -w~o. 



(12) 



(13) 



Presence of these constraints is simply due to the fact that the equations for W and A 
are of the first order in time. The basic equal-time Poisson bracket relations are 

{W(p,t),n w (p',t)} = 8(p-p'), (14) 

{A(p,t),n A (p',t)} = $(p-p')- (15) 



In Dirac's terminology [22], the constraints in Eqs. (12) and (13) are of the second class, 
since {xi (p, t), X 2 (p' > 0/ = ~8(p - p ) which does not vanish weakly. Then, to 
eliminate these second-class constraints, it is conventional to introduce the Dirac 
bracket [22] defined by 



{A(t\ B(t)}* = {A(t), B(t)} 

-j^\\dpdp' {A(0, Xi(P> t)\C u (p, P' ){Xj(p' , t), B(t)}, (16) 

ij=i 
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where A and B are functionals of W and A. In this equation, C ij (p, p' ) are quantities 

satisfying the equations: J^J^' {*,■(/>, 0, , t)}C kj (p" , p' ) = 5 u 5(p- p' ). 

In the present case, C u (p,p') are given as follows: C n (p, p') = C 22 (p, p') = 0, 
C 12 (p, p' ) = -C 21 (p, p') = 8(p- p'). Therefore, we obtain the basic relation 



{W(p,t),A(p',t)Y = 8(p-p'), (17) 



which implies that W and A are canonically conjugate to each other with respect to the 
Dirac bracket. 

The Hamiltonian is given by 



= \dp \K— W-D 



3A„. ^ASW) 



dp d p dp 



and is clearly a constant of motion. Using Eqs. (17) and (18), we can ascertain that time 
evolution of W, that is, 
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^M^O, »).«}'. (19) 
at 



precisely reproduces the Fokker-Planck equation in Eq. (1). The equation for A is also 
described in a similar form. 

Next, we consider the generator of the dilatation transformation 



G=\dp p^W. (20) 
d p 



This quantity gives rise to the following relations: 



{G(t),W(p,t)}* = ^-[pW(p,t)], (21) 
a p 



{G(t),A(p,t)Y=p^P^-. (22) 

dp 



Therefore, the finite transformations are expressed as 



exp {(In A) G (?)}* W (p, t) exp {-(In X) G (t)}* 



= e (lnX) dp " W(p,t) = XW(?ip,t), (23) 
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exp{(lnA)G(0}* Mp, exp{-(lnA)G(0} 



= e 



dp 



A(p, t)=A(Ap, t), 



(24) 



where A is a positive constant and e [A] Be~ [A} = B + {A,B}* + (1 / 2\){A,{A,B}*}* + ■ ■ ■. 
Normalization of W is kept unchanged, whereas the auxiliary field, A, need not be 
normalized. The Dirac bracket relation in Eq. (17) is preserved, as it should be, since the 
transformations are canonical. 

Now, we are at the position to discuss the dilatation symmetry of the system. Such a 
symmetry is characterized by the equation 



This invariance principle may tell us under what conditions the Fokker-Planck equation 
in Eq. (1) admits a scale-free solution. After some calculations using Eqs. (18), (21), and 
(22), we obtain 



{G, H}*=0. 



(25) 



{G, H}* 




-D 



dp{ dp) 



d f dW 



(26) 
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Therefore, invariance under the dilatation transformation gives rise to the condition 



ir tit dK w r, dW d 

KW-P — W-D— + P — 
dp dp dp 







< dw) 






P a 




d p 





= 0. 



(27) 



We note that here W need not be stationary. We rewrite Eq. (27) as 



— f- 
d -lnW= 1 dp ^ P > 



dp p d ( 



d p 



2 

VP J 



(28) 



The asymptotic power-law behavior of W means that 



W(p, t) 



a{t) 



(29) 



holds for large values of p, where a(t) is a positive function of time and the exponent, 
a , is assumed to be a positive constant independent of time. Then, d In W I d p — a I p , 
and Eq. (28) gives the condition 



K(p) , D(p) 



= c, 



(30) 
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where c is a constant. Eq. (30) is our main result. Clearly, Eq. (7) can satisfy Eq. (30) in 
some circumstances including Eq. (2), but the latter is more general than the former. We 
again emphasize that Eq. (30) is free from the assumption of W to be stationary. 

To see the meaning of c in Eq. (30), it is necessary to consider the law of time 
evolution, i.e., the Fokker-Planck equation. Let us take Eq. (1) without assuming Eq. (2). 
Substituting Eq. (29) into Eq. (1) and using Eq. (30), we find 



da(t) , , . ,„ . 

— — ~ c{0 -\)a{t), (31) 

dt 



or its solution 



a(t) ~ a(0) e c{a ~ l)t . (32) 



Thus, we see that c in Eq. (30) is related to the asymptotic factor in Eq. (29) as follows: 



1 , a(t) 

c In . (33) 

o-\ a(0) 



In conclusion, we have developed a general method to assess the asymptotic scale- 
free nature of the solution of the Fokker-Planck equation. We have seen, in the special 
case when the distribution is stationary, how this method can reveal the symmetry 
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principle underlying Lutz's result on the connection between anomalous transport in the 
optical lattice and Tsallis statistics in the high-energy regime. We have also derived the 
condition in the general nonstationary case, under which the solution distribution 
becomes scale-free. 

In the present study, we have treated only the linear Fokker-Planck equation. It has 
been shown [23] that a certain kind of nonlinear generalizations of the Fokker-Planck 
equation, which appear in various physical situations, can also yield solutions that decay 
as a power law. Applications of the present method to such generalized equations and 
systems will be discussed elsewhere. 

This work was supported in part by the Grant-in- Aid for Scientific Research of Japan 
Society for the Promotion of Science. 
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